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We show that the quantum dynamics of interacting and noninteracting quantum particles are 
fundamentally different in the context of solving a particular computational problem. Specifically, we 
consider the graph isomorphism problem, in which one wishes to determine whether two graphs are 
isomorphic (related to each other by a relabeling of the graph vertices), and focus on a class of graphs 
with particularly high symmetry called strongly regular graphs (SRG's). We study the Green's 
functions that characterize the dynamical evolution single-particle and two-particle quantum walks 
on pairs of non-isomorphic SRG's and show that interacting particles can distinguish non-isomorphic 
graphs that noninteracting particles cannot. We obtain the following specific results: (1) We prove 
that quantum walks of two noninteracting particles, Fermions or Bosons, cannot distinguish certain 
pairs of non-isomorphic SRG's. (2) We demonstrate numerically that two interacting Bosons are 
more powerful than single particles and two noninteracting particles, in that quantum walks of 
interacting bosons distinguish all non-isomorphic pairs of SRGs that we examined. By utilizing high- 
throughput computing to perform over 500 million direct comparisons between evolution operators, 
we checked all tabulated pairs of non-isomorphic SRGs, including graphs with up to 64 vertices. (3) 
By performing a short-time expansion of the evolution operator, we derive distinguishing operators 
that provide analytic insight into the power of the interacting two-particle quantum walk. 

PACS numbers: 03.67.Lx,05.40.Fb,02.10.Ox,03.67.Ac 



I. INTRODUCTION 

Random walks have been applied successfully to many 
problems in physics, as well as in many other disciplines, 
stretching from biology to economics P~H5]. Most of 
these applications use classical random walks (CRWs), 
in which quantum mechanical principles are not consid- 
ered. However, more recently, researchers have found 
that CRWs and quantum random walks (QRWs) can 
exhibit qualitatively different properties [6H1]. From a 
standpoint of algorithms research, these disparities lead 
to situations in which algorithms implemented by QRWs 
can be proven to run faster than the fastest possible clas- 
sical algorithm [TJHTFj]. 

Besides being useful as theoretical models, simple 
QRWs have already been experimentally implemented 
in externally driven cavities [16J . arrays of optical traps 
OH HE], NMR systems QU, and ion traps [20l EIJ . This, 
coupled with new ideas for realistic physical implemen- 
tations of non-trivial walks [22J, indicates that studying 
algorithms cast as QRWs might lead to algorithms that 
are both powerful and experimentally viable. 

Although QRWs are universal and therefore in princi- 
ple can be used to implement any quantum algorithm 
[23) . in practice some information-theoretic problems 
lend themselves to a QRW approach more easily than 
others. Many of the problems that have been investigated 
are expressed naturally using graphs, sets of vertices and 
edges, with each edge connecting two vertices. For exam- 
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pie, the vertices of a graph might be taken to represent 
individuals, and the edges between them might indicate 
friendship. A question one could ask is whether there is 
a subset of friends that are isolated from the rest of the 
group, which translates to the graph-theoretic question 
of whether the graph is disconnected. Once the question 
to be answered is posed as a question about a graph, it 
is investigated by constructing a quantum Hamiltonian 
from that graph. The dynamics of the system is then 
analyzed using quantum mechanics, and is used to make 
statements about the original graph, hopefully giving in- 
sight to the answer of the original problem. 

This paper addresses the graph isomorphism (GI) 
problem, where, given two graphs, one must determine 
whether or not they are isomorphic (two graphs are iso- 
morphic if one can be obtained from the other by a re- 
labeling of the vertices). Although many special cases of 
GI have been shown to solvable in a time that scales as 
a polynomial of the number of vertices, the best general 

classical algorithm to date runs in time O (c N 1 ^s-^V 

where c is a constant and N is the number of vertices in 
the graphs being compared [2"3] . 

GI has several properties analogous to those of factor- 
ing. First, though it appears to be difficult, it is felt that 
it is unlikely to be NP-complete, since otherwise many 
complexity classes believed to be distinct would collapse 
|43j . Second, both GI and factoring can be viewed as 
hidden subgroup problems. In GI, one is looking for 
a hidden subgroup of the permutation group, while in 
factoring, one is looking for a hidden subgroup of the 
Abelian group. The success of Shor's polynomial-time 
algorithm for factoring [25] has led several groups to in- 
vestigate a hidden-subgroup approach to GI. However, 
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the obstacles facing such an approach have been shown 
to be formidable [26| [27] . 

Researchers have also recently attacked GI using var- 
ious methods inspired by physical systems. Rudolph 
mapped the GI problem onto a system of hard-core 
atoms [28 . One atom was used per vertex, and atoms 
i and j interacted if vertices i and j were connected by 
edges. He showed that pairs of non- isomorphic graphs 
exist whose original adjacency matrices have the same 
eigenvalues, while the induced adjacency matrices de- 
scribing transitions between three-particle states have 
different eigenvalues. Gudkov and Nussinov proposed 
a physically-motivated classical algorithm to distinguish 
non- isomorphic graphs [29] , Shiau et al. proved that 
the simplest classical algorithm fails to distinguish some 
pairs of non-isomorphic graphs |30| and also proved that 
continuous-time one-particle QRWs cannot distinguish 
some non-isomorphic graphs [30]. Douglas and Wang 
modified a single-particle QRW by adding phase inho- 
mogeneities, altering the evolution as the particle walked 
through the graph [31 j - Their approach was power- 
ful enough to successfully distinguish many families of 
graphs considered to be difficult to distinguish, including 
all families of strongly regular graphs they tried. Most 
recently, Emms et al. used discrete-time QRWs to build 
potential graph invariants |321 133) . Through numerical 
spectral analysis, they found that these invariants could 
be used to distinguish many types of graphs by breaking 
the eigenvalue degeneracies of many families of graphs 
that arc difficult to distinguish. 

In addition to studying single-particle QRWs, Shiau 
et al. [30] performed numerical investigations of two- 
particle QRWs and presented evidence that interacting 
Bosons can distinguish non-isomorphic pairs that single- 
particle walks cannot. There, it was also found numeri- 
cally that two-Boson QRWs with noninteracting particles 
do not distinguish some non-isomorphic pairs of graphs. 
In contrast to the approaches in 31-33], the two-particle 
QRW algorithm does not lower the symmetry of the sys- 
tem. 

In this paper, we extend the results in [30] on two- 
particle quantum walks in several ways. First, we prove 
analytically that quantum walks of two non-interacting 
Bosons always fail to distinguish non-isomorphic pairs of 
strongly regular graphs (SRGs). This result is surprising, 
since it has been shown in [34-35] that non-interacting 
Boson QRWs on graphs can give rise to effective statis- 
tical interactions, which significantly alter the dynamics 
of the system. Second, we show that analysis of the be- 
havior of non-interacting Fermions requires a more sub- 
tle treatment than was done in [3D]; the result in [30] 
that some non-isomorphic SRGs could be distinguished 
by two noninteracting Fermions and not by two noninter- 
acting Bosons arose because of an ambiguity in the choice 
of basis. When the ambiguities involved with the basis 
choice for Fermions are accounted for, non-interacting 
Fermions fail to have any advantage over non-interacting 
Bosons. Third, we expand on the initial numerical results 



in [30] , exhaustively verifying, where only sampling was 
used before, that two-particle interacting Boson walks 
distinguish all the non-isomorphic pairs of SRGs with up 
to 64 vertices that have been tabulated. To accomplish 
this, we used high-throughput computing techniques to 
perform over 500 million comparisons between evolution 
operators of graphs. Finally, we examine the small-time 
expansion of the evolution operator, and use the two- 
particle interacting evolution to derive candidates for 
graph invariants, which appear in the fourth and sixth 
orders in time. 

Our results demonstrate unambiguously that two- 
particle Bosonic quantum walks have more computa- 
tional power if the particles are interacting, because inter- 
acting walks can be used to distinguish non-isomorphic 
graphs that noninteracting particles cannot. 

The paper is organized as follows. Section|n]introduces 
relevant background and definitions to QRWs on graphs, 
including a brief overview of the strongly regular graphs 
(SRGs) on which the algorithms are tested and also a 
review of the one-particle QRW algorithm considered in 
[50] . Section III proves that QRWs of two non-interacting 



Bosons do not distinguish non-isomorphic SRGs with the 



same family parameters. Section IV analyzes the QRW 
of two non-interacting Fermions. It shows that improper 
basis choice can lead to false-distinguishing, and that 
when the basis is chosen consistently, QRWs with two 
noninteracting fermions are unable to distinguish some 
pairs of non-isomorphic SRGs. Section [V] shows through 
exhaustive simulation that all tabulated families of SRGs 
are successfully distinguished by a two hard-core Boson 



QRW. In section VI a short-time expansion of the evolu- 
tion operator is computed, and distinguishing operators 
present in the two hard-core Boson QRW are identified. 
Finally, Section |VII| summarizes and discusses the pos- 
sible implications of our results for the development of 
algorithms based on interacting QRWs for distinguishing 
non-isomorphic graphs. 



II. BACKGROUND & DEFINITIONS 



This section describes the background and definitions 
necessary to study QRWs on graphs. First, we introduce 
the graph-theoretic concepts we will need, including the 
notions of of the adjacency matrix and the spectrum of a 
graph. Then, we consider how to use these tools to con- 
struct a physical process through the definition of Hamil- 
tonian and Green's functions, for both one and two par- 
ticles. Next, we detail the relevant properties of SRGs. 
Finally, we review the method used in 30J to use the 
properties of SRGs to show that the Green's functions of 
single-particle QRWs do not distinguish non-isomorphic 
SRGs with the same family parameters. 
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A. Constructing walks on graphs 

In this section, we describe how to form QRWs on 
graphs. A graph G — (V, 25) is a set vertices V and 
edges E. The vertices are usually labeled by inte- 
ger indices, and the edges are unordered pairs of ver- 
tices. Two vertices that share an edge are called con- 
nected, while two vertices that do not are called dis- 
connected. The total number of neighboring vertices of 
a particular vertex is called its degree. For example, 
the graph G = ({1,2,3,4}, {(1,2), (2,3), (3,4), (4,1)}) 
is a cycle graph with four vertices. Two graphs are 
isomorphic if they can be made identical by relabel- 
ing their vertices. For example, the graph H = 
({1,2, 3, 4}, {(1,3), (3, 2), (2, 4), (4,1)}) is isomorphic to 
G, since after relabeling 2 o 3, the two graphs are equiv- 
alent. 

Graphs are conveniently expressed algebraically as ad- 
jacency matrices. An adjacency matrix A of a graph with 
N vertices is an N x N matrix in the basis of vertex la- 
bels, with Aij = 1 if vertices i and j are connected by an 
edge, and zero otherwise. The adjacency matrix for the 
graph G is 




(1) 



The spectrum of a graph is the eigenvalue spectrum of its 
adjacency matrix. The spectrum of G is {—2, 0, 2}, with 
two-fold degenerate. 

To form a QRW on a graph, we first define a Hamil- 
tonian. We will use the Hubbard model, with each site 
corresponding to a graph vertex. A particle can make a 



transition between two sites if the associated vertices are 
connected. In addition, if two walkers happen to simul- 
taneously occupy a site, we impose a double-occupation 
energy cost U . Our Hamiltonian is 



H 



J2 A 



u 

~2 



ejci 



c\c t 



(2) 



where c and c< are Boson or (spinless) Fermion creation 
and annihilation operators. If we restrict ourselves to 
single-particle states, we find matrix elements 



(*|H|j) 



-A, 



(3) 



Hence, we can easily identify the a single-particle Hamil- 
tonian 



H 1P = A. 



(4) 



Similarly, we can define two-particle Hamiltonians by 
their matrix elements. For these, we need to use either 
Bosonic or Fermionic basis states. The Boson states are 



i = j 



and the Fermion states are 



1 

71 



m-m)- 



(•5) 



(6) 



We now restrict ourselves to two particles, where we de- 
fine H2B and H2F to be the two-particle Boson and 
Fermion Hamiltonians. These are special cases of eq. 
[2] with matrix elements 



-SikAji - SjiA ik - SuAjk - S jk Au : j and k ^ I 
} {ij\H 2B \kl) B = B <y|H \kl) B = { US lk : i = j and k = I , 

(SikAji + S^Aik + S u Ajk + d jk A u ) : i = j xor k = I 

F (ij|H 2 _F \kl) F = F {ij\H\kl) F = A ik Sji + AjiSik - AuSjk - AjkSu : i ^ j and k ^ I, (7) 



respectively, where the matrix elements are found di- 
rectly from Eq. [2] through the application of appropriate 
commutation relations. 

From each of these Hamiltonians, we define the QRW 
time-evolution operator as 



U = e 



-iiH 



(8) 



where we set h = 1 for notational convenience. 

To study the dynamics of the system, we define the 
two-particle Green's function (GF), which relates the 
wavefunctions at a time t to those at time t = 0. For 



two-particle position states \ip) and \ip'}, it is 

g (v(o)>'W) = (m i/(t)) = (v>(0)| u |v'(o)> . (9) 



Letting (V>,VO run over a complete two-particle basis, 
0(^>(O) , ip'(t)) considered at a fixed time provides us with 
a set of N 2 (N + l) 2 /4 complex numbers for Bosons or 
N 2 (N — l) 2 /4 complex numbers for Fermions. These 
lists completely characterizes the dynamics of the sys- 
tem. Hence, when we analyze QRWs, we say that a par- 
ticular scheme distinguishes two graphs if their Green's 
functions supply us with two different lists. 
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B. Strongly Regular Graphs 

Our major results in this paper focus on algorithms 
over strongly regular graphs (SRGs), which are difficult 
to distinguish. In this section, we develop the properties 
of SRGs we will need for our later analysis. 

A SRG is a graph in which (a) all vertices have the 
same degree, (b) each pair of neighboring vertices has the 
same number of shared neighbors, and (c) each pair of 
non-neighboring vertices has the same number of shared 
neighbors. This definition permits SRGs to be cate- 
gorized into families by four integers (N, k, A,/x), each 
of which might contain many non-isomorphic members. 
Here, N is the number of vertices in each graph, k is 
the degree of each vertex (k-regularity), A is the number 
of of common neighbors shared by each pair of adjacent 
vertices, and [i is the number common neighbors shared 
by each pair of non-adjacent vertices. 

Using the stringent constraints placed on SRGs, one 
can show that, regardless of size, the spectrum of any 
SRG only has three distinct values |37j : 

A = -k, (10) 

which is non-degenerate, and 

Xi, 2 = -1(\-»±Vn), (11) 

which are both highly degenerate. Both the value and 
degeneracy of these eigenvalues depend only on the fam- 
ily parameters, so within a particular SRG family, all 
graphs are cospectral [37]. Further, the spectra of the 
two-particle Hamiltonians formed from SRG adjacency 
matrices, as described in Eq. [7J are also highly degen- 
erate. As shown in figure [j] for the family (16,6,2,2), 
the interacting case gives us the largest number of dis- 
tinct energy levels. These highly degenerate spectra are 
one reason why distinguishing non-isomorphic SRGs is 
difficult — it is known that distinguishing non-isomorphic 
graphs with spectra with bounded degeneracy can be 
done with polynomially bounded resources |38j . 

The adjacency matrix of any SRG satisfies the useful 
relation [3"Tj : 

A 2 = (k - fj)I + fjj + (A - fi)A, (12) 

where I is the identity and J is the matrix of all ones 
(J y = 1 for all Since J 2 = iVJ and I 2 = I, this 

forms a three-dimensional algebra, and we can write 

A" = a n I + (3 n J + 7n A, (13) 

where a, /3, and 7 are functions only of the family pa- 
rameters. That is, all SRGs of the same family have the 
same coefficients. 

Although many SRGs are known there are sub- 
stantially fewer tabulated families with more than one 
non-isomorphic member. Complete and partial families 
of SRGs have been tabulated through combinatorial tech- 
niques [101 HI], which we use in sec. |V|to perform nu- 
merical tests of our algorithms. 
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FIG. 1: The energy spectra for several types of QRWs on 
the SRG family (16,6,2,2), where the parameters (N , fc, A, /z) 
define an individual SRG family. N is the number of vertices 
in each graph, k is the degree of each vertex (k-regularity), 
A is the number of of common neighbors shared by each pair 
of adjacent vertices, and fi is the number common neighbors 
shared by each pair of non-adjacent vertices. Both graphs 
in the family have the same spectra. In all four panels, the 
Hamiltonian used is the Hubbard model (eq. [2j) , with U = 
in the noninteracting cases and U — > 00 for the hard-core 
Bosons. The degeneracies are given to the right of each level. 



C. Review of one-particle algorithm 

In this section, we review the method used in Ref. [30] 
to prove that a single-particle QRW cannot distinguish 
non-isomorphic members of SRG families. Formally, this 
means that we must show that any two SRGs of the same 
family parameters have the same single-particle Green's 
functions Q\p. 

First, we consider the adjacency matrix A, and sup- 
pose it belongs to the strongly regular graph family 
(N,k,X,fi). Then, by Eq. gj we know that the single- 
particle hamiltonian is Hip = —A, and so by Eq. [U the 
QRW evolution operator is Uip = e ltA . But since A is 
a SRG, we make use of the algebra defined in Eq. [T3]to 
write 

Uip = aI + /3J + 7 A, (14) 

where the coefficients depend only on the family param- 
eters. 

Following Shiau et al., we investigate the relevant 
Green's functions, Gip(i,j) = (i\TJip\j). We first con- 
sider the diagonal elements, each of which contains a con- 
tribution of a from I and f3 from J. Note that there is no 
contribution from A, because it is entirely off-diagonal. 
Hence, the N diagonal Green's functions are all equal 
to a + p. For the off-diagonal elements, I never con- 
tributes and J always does. However, A contributes to 



5 



only some of the elements. More precisely, each column 
of A contains exactly k ones (entirely in the off-diagonal), 
since each vertex is of degree k. Hence, there are kN off- 
diagonal Green's functions of the form f3 + 7, and the 
remaining iV 2 — N — kN are equal to /3. 

As can be seen from Eq. [13J a, /3, and 7 all de- 
pend only on the family parameters. Therefore, the one- 
particle evolution for any graph in the same family will 
have the same GFs, and the algorithm based on single- 
particle quantum evolution fails to distinguish any non- 
isomorphic SRGs that are in the same family. 



III. PROOF THAT QRWS WITH TWO 
NON-INTERACTING BOSONS DO NOT 
DISTINGUISH NON-ISOMORPHIC SRGS IN 
THE SAME FAMILY 

We now develop exact expressions for the time evolu- 
tion operators for two non-interacting Bosons, and sub- 
sequently show that this evolution cannot be used to dis- 
tinguish non-isomorphic SRGs in the same family. Al- 
though this result may seem expected, recent efforts [531 - 
155] have demonstrated that non-interacting QRWs on 
non-translationally invariant graphs lead to effective, sta- 
tistical interactions, resulting in rich physical phenomena 
such as Bose-Einstein condensation. The method used 
here is analogous to that used in Ref. [3D] to show that 
one-particle QRWs cannot distinguish non-isomorphic 
SRGs from the same family, but with a more complex 
implementation. 

First, we note that we may write the Hamiltonian for 
any two-Boson QRW as 

H 2B = -i(I + S)(A0A) + [/R, (15) 

where A©A = A(S>I + I<X>Aisa Kronecker sum, the 
matrix special case of a direct sum, and 

s = E ly') 0*1 . r = E l <4 > <"l • ( 16 ) 

The demonstration that Eq. [15] is equivalent to Eq. [7] is 
given in appendix [A] For noninteracting Bosons [7 = 0, 
and the evolution operator is 



= E^ K (I+s)(AeA) ' (17) 

71 = ' ^ ' 

where U 2 b is shorthand to refer only to the non- 
interacting case. Now, by the definitions of the Kronecker 
sum and S, it is easy to show that [(I + S) , (A © A)] = 0. 
Hence, 

00 , . n 

U2B = E^K (I + S)"(A©Af. (18) 

n=0 ' ^ ' 



But note that S 2 = I, so (I + S)" = 2™- 1 (I + S). It 
follows that 

U 2S = i(I + S)e ltA0ltA . (19) 

Since each matrix A in the Kronecker (direct) sum 

exponentiates in its own product space, we can write 
e A0B _ g A g B ^ wn j cn i eac i s to 

U 2B = - (I + S) e UA ® e ltA 

= i(I + S)U lp ®U lp , (20) 

where the single-particle evolution operator Uip is de- 
fined in Eq. [8] for H = —A. Since the Boson states are 
symmetric under particle interchange, we have matrix 
elements 

B (ij\V 2B \kl) B = B (ij\V 1P <g> Ulp \kl) B . (21) 
Expanding this using eq. |14[ we have 

B (ij\U 2B \kl) B = B (ij\(a 2 I(Z>I + [3 2 J ® J 

+ 7 2 A ® A + a/3 (J J) + 07 (A © A) 
+ #y (J ® A + A ® J) ) \kl) B . (22) 

Now that we have determined the matrix elements of 
U 2 b , we can work out all the cases for eq. [22j which are 
the GFs of the system. We find that there are 22 pos- 
sible values, each of which can be written as an explicit 
function of a, ft, and 7. Since the values of the matrix el- 
ements are all only functions of SRG family parameters, 
they are the same for all graphs in the same SRG family. 
To show that the GFs are the same across a SRG fam- 
ily, we also need to show that the number of occurrences 
of each value is also a function only of SRG parameters. 
In Appendix [B] we count all the types of these GFs in 
terms of SRG family parameters, with the results shown 
111 table |H 

Because we have shown that both the values and num- 
ber of occurrences of all of the two-particle Green's func- 
tions for noninteracting Bosons can be written in terms of 
the family parameters N, k, A, and fi, we demonstrated 
that two noninteracting Bosons cannot distinguish non- 
isomorphic SRGs of the same family. 

IV. ANALYSIS OF THE NONINTERACTING 
TWO-FERMION EVOLUTION FOR SRGS 

In this subsection we consider the analogous evolution 
generated by two non-interacting Fermions. This analy- 
sis is more complicated than for Bosons because changing 
the two-particle basis can introduce sign changes. If this 
sign ambiguity is not accounted for properly, the algo- 
rithm may falsely distinguish two graphs that are actu- 
ally isomorphic. 
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TABLE I: This table enumerates the Green's functions for the QRW with two non-interacting Bosons on a SRG with family 
parameters (TV, fc, fi, A). The Hamiltonian considered is H_b = — 1/2 ■ (I + S) (A © A), where S swaps the two particles and A is 
the adjacency matrix of the graph. Because the A, I, and J form an algebra, the parameters a, 0, and 7 are the same for every 



graph in an SRG family (Eq. (13 1. The evolution operator for noninteracting Bosons, Us = 1/2 (I + S) U 



U 1 



contains 



terms bilinear in I, J, and A, with coefficients that can be written in terms of a, 0, and 7. The GFs, formed by taking matrix 
elements of Ub (Eq. |9J, are divided into symmetry classes (a, b), where a indicates the number of distinct basis indices and 6 
is the number of indices shared between the left and right sides. For example, (34| B Ub |24)„, where B indicates identical 
Bosons on vertices i and j, falls into the symmetry class (3, 1), since it has three distinct indices (2,3,4) and the left and right 
side have one index in common (2). Since the total number of entries with any given particular value can be counted in terms 
of numbers that are constant for a given set of family parameters, the GFs generated by non-isomorphic members of the same 
SRG family must have the same values and the same degeneracies. 



Element Class 


Value of Element 


Number of Occurrences 


(4,0) 


40 7 + 2 7 2 + 2/3 2 


1/4 ■ TV (k 2 (fi + 1) + k (A 2 - X(n + 2) + fi - 1) - 2(N - 




3/37 + 7 2 + 2/3 2 


Nfi(N - k-l)(k + X- n) 




20 7 + 20 2 + 7 2 


l/(2fc) • [N(N — k — 1) (k 3 — 2k 2 fi + (N- 1)/j, 2 )] 




20 7 + 2/3 2 


1/k ■ [N(N - k - 1) (fc 3 - fe 2 (2/i + 1) + (TV — l)/i 2 )] 




07 + 20 2 


1/k ■ [N(k - N + l)(fe - fi)(k(2k — TV + 2) — N/i + fj,)] 




20 2 + 7 2 


l/(4fc) • [TV(TV - fc - l)(fc(-3fcTV + fc(3fc + 8) 
+TV 2 — 5TV + 6) — 2fc(fc + l)n + (TV - l)^ 2 )] 




Subtotal: 


1/4 • TV(TV - 1)(TV - 2)(TV - 3) 


(3,0) 


i/y/2 ■ 07 + 2/v/2 • 2 + 2/V2 ■ 7 2 


kN(k - A - 1) + fcTVA 




2/^2 • /9 7 + 2/^2 ■ 2 


2fcTV(TV - fc - 1) 




2/V2- 2 


TV(fc-TV + l)(fc-TV + 2) 




Subtotal: 


TV(TV- l)(TV-2) 


(3,1) 


a0 + a 7 + 3/37 + 20 2 + 7 2 


fcTVA 




a/3 + 2/3 7 + 2/3 2 + 7 2 


TV(TV - 1 - k)n 




a/3 + 07 + 207 + 20 2 


2N(N- l-k)fi 




a0 + a 7 + 7 + 20 2 


fcTV(-2fc + TV + A) 




a0 + 07 + 20 2 


2fcTV(— 2fc + TV + A) 




a0 + 20 2 


TV(1 + fc - TV)(2 + 2k - TV - fj) 




Subtotal: 2 


TV(TV- l)(TV-2) 


(2,0) 


207 + 2 + 7 2 


fcTV 




2 


TV(TV-fc- 1) 




Subtotal: 


TV(TV- 1) 


(2,1) 


2/V2- (q0 + a 7 + 07 + 2 ) 


2fcTV 




2/V2- (a0 + 2 ) 


2TV(TV - fc - 1) 




Subtotal: 


2TV(TV- 1) 


(2,2) 


a 2 + 2a0 + 207 + 20 2 + 7 2 


1/2- fcTV 




a 2 + 2a0 + 20 2 


1/2 ■ TV(TV - fc - 1) 




Subtotal: 


l/2-TV(TV-l) 


(1,2) 


a 2 + 2a0 + 2 


N 


Total: 




1/4- TV 2 (TV + l) 2 



The Hamiltonian for the two-Fermion QRW, 'H.2F^ is logic we did for the Bosons but let t —> —t, 



H 2F = - (I — S) (A © A) 



(23) 



V 2F (-t) = U 2F 



(24) 



where, again, / is the identity, S is the operator that 

swaps the two particles, A is the adjacency matrix of the This way, the single-particle evolution operator Uip still 
graph, and © denotes a direct sum. We follow the same has time running forwards. The matrix elements of U2_f 
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TABLE II: Enumeration of the Green's functions for QRWs of two non-interacting Fermions. The Hamiltonian for two 
noninteracting Fermions is H2F = 1/2 • (I — S) (A© A), where S swaps the two particles and A is the adjacency matrix 



of the graph. 
U F = 1/2 (I- 



Hence, the two-Fermion evolution operator, defined with t 



-t to keep U running with forward time, is 



S)U X 



U .It contains terms bilinear in I, J, and A, with coefficients written as combinations of a, /3, and 



7, where the parameters a, /3, and 7, are defined in Eq. (141. We divide the matrix elements {ij\ F U2F \kl) F of the GFs into 



classes (a, b), where a indicates the number of distinct indices and b is the number of indices shared between the left and right 
sides. For example, (34j F Uf |24} f falls into class (3, 1), since it has the three distinct indices (2,3,4) and the left and right side 
have the index (2) in common. The ± next to some of the element values indicates that the count applies to all elements with 
the given magnitude; this is done because the number of elements with each sign depends on the choice of two-particle basis. 
The total number of matrix elements with a given absolute value can be expressed in terms of the SRG family parameters. 
Therefore, we conclude that the matrix elements of the GFs of two non-isomorphic members of the same SRG family must be 
equivalent up to sign differences. 



Element Class 


Value of Element 


Number of Occurrences 


(4,0) 





l/(4fe) • n( - 6k 4 + 2fc a (5n + 6/z - 7) - 4fc 2 (n - l)(n + 3/i - 2) 
+fc(n - 1) ((n - 5)n - 6/i 2 + 6) + 6(n - 1) V) 




±(/? 7 + 7 2 ) 


Nfj, ■ (N - k - l)(k + A - fx) 




± ( 7 2 + 2/3 7 ) 


1/k ■ [N(N - k - 1) (fc 3 - 2fe 2 /x +(N - l)n 2 )] 




±Pl 


1/k ■ [N(k - N + l)(k - n){k{2k - N + 2) - Nfi + fi)] 




Subtotal: 


1/4 ■ N(N - i)(N - 2)(N - 3) 


(3,1) 


± (a/3 + a?7 — /?7 — 7 2 ) 


fciVA 




± (a/3 - 2/3 7 - 7 2 ) 


N(N - 1 - k)n 




± (a/3 + 07) 


2N(N- 1 - k)fj, 




± (a/3 + «7 + /3 7 ) 


kN(-2k + N + A) 




± (a/? - M) 


2kN(~2k + N + A) 




±a/3 


iV(l + k — N)(2 + 2k — N — n) 




Subtotal: 


N{N - 1)(N -2) 


(2,2) 


a 2 + 2a/3 - 2^7 - 7 2 


1/2 -kN 




a 2 + 2a/3 


l/2-N{N-k-l) 




Subtotal: 


1/2 -N(N- 1) 


Total: 




1/4 ■ N 2 (N — l) 2 



are given by 

F {ij\V 2 F \kl) F = p{ij \ [a 2 l ® I + /3 2 J ® J + 7 2 A ® A 

+ a/3 (J 8 J) + a-f (A ® A) 

+ /3 7 (J® A + A® J)) \kl) F . (25) 

We now show that there are sign ambiguities in Uf 
that arise from the choice of basis states that one uses 
when converting a graph to an adjacency matrix. As an 
example, consider the two isomorphic graphs shown in 
figure [2j which have adjacency matrices 



(26) 



We wish to write down the two-particle Hamiltonians 
using eq. [7] However, we must first pick a basis. That 
is, for each pair of sites \ij) = \ji), we are free to pick 
either ordering, but we cannot choose both. Supposing 





1 








1 




i 







■ - 







i 




1 
















Case A: 



Case B: 



FIG. 2: Two clearly isomorphic graphs. Graph A differs from 
graph B only by the labeling of vertices 2 and 3. Despite 
this, some matrix elements of two-particle Fermion evolution 
operators Uf of the two graphs have opposite signs. This 
result implies that using the numerical values of these matrix 
elements produces a false-positive: two isomorphic graphs are 
falsely distinguished. 



we pick {|12) , |13) , |23)}, we get 





(0 


1 






( 1 




Ha = 


1 





!) 


, Hg = 


1 


:) 




\° 


1 






v -i 





(27) 



Forming the evolution operators U 



itH , we have 



/ / cc^v^t) j\ isin(V2t) I cosf^t) i \ \ 



U„ = 



and 



i(V2t) 



U, 



V2 
cos(\/2t) 



,(V2t) 

dn(\/2t) 



COS i\Fit 



V2 

isin(V2t) / cosfv^t) 



Fit) 



V2 



(28) 



V 



V2 



/2t) 




The values of the matrix elements of these two evolution 
operators clearly have sign differences. If we had taken 
into account that the second and third labels had been 
switched on graph B, we would have chosen the basis to 
be {|12) , 1 13} , |32)}, fixing the factors of negative one. 
Unfortunately, in a situation where we are handed two 
graphs and asked whether or not they are isomorphic, 
we do not know a priori what the correct basis choice 
should be for proper testing. Moreover, because the num- 
ber of possible basis choices is 2 N , checking all of them 
is not feasible. We adopt here a simple strategy that 
eliminates this dependence of the sign on the choice of 
basis that arises for more than one fermion; we compare 
the absolute value of the elements, rather than the ele- 
ments themselves. The absolute values of all the elements 
and their degeneracies are shown in Table [TT] The abso- 
lute values and degeneracies can be expressed as explicit 
functions of family parameters, so we conclude that two 
non-interacting Fermions, as well as two non-interacting 
Bosons, fail to distinguish non-isomorphic SRGs from the 
same family. The enumeration of the three classes of ma- 
trix elements allowed by Uf are listed in table [TTJ where 
the ± symbol indicates that the count given is the total 
of elements of either sign. 



V. NUMERICAL TESTING OF EVOLUTIONS 
OF RANDOM WALKS OF INTERACTING 
PARTICLES 

In the preceding two subsections we have proven that 
QRWs with two noninteracting particles are not useful 
for distinguishing non-isomorphic SRGs from the same 
family. In this section we perform a systematic investi- 
gation of the ability of QRWs of two interacting Bosons 
to distinguish non-isomorphic SRGs. We go beyond the 
sampling performed in [30 by exhaustively checking the 
two-Boson interacting QRW on all tabulated SRG fami- 
lies with more than one member. Our work shows that 
this walk succeeded in all trials preformed, including suc- 
cessfully distinguishing the (36,15,6,6) SRG family, which 



TABLE III: Numerical simulations of QRWs of two hard- 
core Bosons on many SRG families with multiple non- 
isomorphic members. The Hamiltonian used was Hs = 
-1/2 ■ (I + S) (A © A) + UK, where S swaps the two par- 
ticles, A is the adjacency matrix of the graph, and R counts 
double-occupation. To evaluate the hard-core limit, we took 
U — > oo. All non-isomorphic graphs in the families indicated 
were compared pairwise, with A is a measure of how differ- 
ent the matrix elements of the evolution operator are, defined 
precisely in eq. |30| When A = 0, the list of matrix elements 
of the two evolution operators being compared have the same 
magnitudes. The minimum values of A were non-zero for all 
pairs of non-isomorphic graphs that were examined. 



SRG family (TV, k, fi, A) 


non-isomorphic members 


minimum A 


(16,6,2,0) 


2 


94.273 


(16,9,4,6) 


2 


2.723 


(25,12,5,6) 


15 


3.636 


(26,10,3,4) 


10 


7.356 


(28,12,6,4) 


4 


27.607 


(29,14,6,7) 


41 


4.017 


(35,18,9,9) 


227 


5.243 


(36,14,4,6) 


180 


2.621 


(36,15,6,6) 


32,548 


1.512 


(37,18,8,9) 


6,760 


4.310 


(40,12,2,4) 


28 


3.065 


(45,12,3,3) 


78 


5.868 


(64,18,2,6) 


167 


2.574 



has 32,548 non-isomorphic members. We used the follow- 
ing procedure for each pair of graphs in each family: 

1. Begin with the evolved (complex) evolution matrix 
U A . 

2. Take the magnitude of each element. 

3. Write all the (real) entries to a list, X^. 

4. Sort the list. 

5. Compare the list using 

A = ^|X A L7']-X s [j]|. (30) 

3 

Note that for any two isomorphic graphs, A = 0, so if 
A / 0, we can conclude that the graphs are not iso- 
morphic. Table [TTT] presents our results, which show that 
QRWs of two hard-core Bosons successfully distinguished 
all pairs of non-isomorphic graphs in all SRG families 
tested. 

The process of numerically checking that each pair of 
non-isomorphic graphs was indeed distinguished by our 
algorithm was computationally intensive. First, in order 



to calculate U 



-itH 



one must diagonalizc H. For the 



N = 64 cases we considered, since the two-particle space 
has dimension 64- (64 + l)/2 = 2080, we needed to diago- 
nalize large, non-sparse matrices. Further, we needed to 
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perform many comparisons to generate all the candidates 
for the minimum A. For example, for the (36,15,6,6) fam- 
ily, one needs to perform 32, 548-32, 547/2 = 529, 669, 878 
comparisons to check each pairwise A. To accomplish 
this, we used high-throughput computing environment 
Condor running on the University of Wisconsin's Center 
for High Throughput Computing cluster. The numerical 
error on A was between 10 -14 and 10 -9 for all families 
we analyzed. 

In addition to our calculations with hard-core Bosons, 
we also investigated non-interacting Fermions numeri- 
cally using the numerical procedure described above. As 
discussed above, comparing absolute values of matrix el- 
ements eliminates the sign discrepancy brought on by 
basis choice. As one expects given the results in section 



IV the result A = is obtained for all cases tested (the 



first six SRC families appearing in table III I . 



VI. 



SMALL-TIME EXPANSION AND POSSIBLE 
DISTINGUISHING OPERATORS 



In this section we attempt to gain more insight into 
the behavior of the QRWs of two interacting Bosons by 
expanding the evolution operator for short times. By 
forming such an expansion and listing all the forms that 
appear, we can investigate which of these operators con- 
tribute to the distinguishing power of the evolution op- 
erator. The operators that contribute are called dis- 
tinguishing operators, and below we work out the first 
one, which appears in the fourth order in time. Wc then 
briefly investigate a sixth-order term, which succeeds in 
some instances when the fourth-order term fails. 

We begin with the exact evolution operator for the 
interacting Boson case, which is 



U = 



A(I+S)(A0A)+[/R) 



(31) 



Expanding as a power series in t, we have 



n=0 



i(I + S)(A0A) + [/R j . c-52) 



Since we know from simulation that this evolution op- 
erator distinguishes SRGs, we expect that there will be 
an order in t at which there are terms that are not func- 
tions of only the SRC family parameters. Numerically, 
as shown in figure [3j we calculate that in the case shown 
these terms first appear at fourth order in time, so we 
endeavor to analyze the matrix elements of the first five 
terms of the expansion for the evolution operator, 

t 2 it 3 t A 

U ~ U - itUi - -U 2 + T U 3 + — . (33) 

We first expand the terms to simplify the evolution op- 
erator, taking A © A = B. 



Ui 



1 




23456789 
Order of time in expansion 



10 11 12 13 



FIG. 3: Small-time expansion comparison for the two non- 
isomorphic SRGs in the family (16,6,0,2) using the interacting 
two-Boson Hamiltonian H 2 b = -1/2 • (I + S) (A © A) + UR 
with U = 50, evaluated with t = 0.01. The evolution op- 
erator was expanded to the different orders in t, and A, a 
measure of differences in evolution operator matrix elements, 
is plotted versus the order of the expansion. The actual 
value of A obtained numerically using the full evolution oper- 
ator is given by the dashed line. The distinguishing operator 
U (A © A) R (A © A) 2 + U (A © A) 2 R (A © A), where A is 
the adjacency matrix of the graph, and R counts double- 
occupation, causes the two graphs to be distinguished at 
fourth order in time. 



U 2 = - (I + S) (f/ 2 R + B 2 - f/RB - LTBR) , (35) 



U 3 = -(I+S) (C/ 3 R + [/B 2 R-(/ 2 BR 

- ?7 2 RB B 3 + [/RB 2 + [/BRB), (36) 

where we used the fact that any term containing the 
product RBR vanishes due to the construction of B. 
Finally, 

U 4 = 1 (I + S) (J7 4 R + J7 2 B 2 R - <7 3 BR - <7B 3 R 
+ C/ 2 RB 2 R - C/ 3 RB C/B 2 RB + 77 2 BRB 



+ L7 2 RB 2 + B 4 17RB 3 17BRB 2 



(37) 



(I + S) (UK - B) , 



(34) 



Numerically, we determined that the fourth order oper- 
ators that successfully distinguished (16,6,0,2), and thus 
could not be counted in terms of SRC family parameters, 
was the combination of ITBRB 2 + £/B 2 RB. If either of 
these operators is removed from the fourth order term, A 
drops to zero. Upon further investigation, graphs with 
more vertices, starting with (25,12,5,6), are not neces- 
sarily distinguished by the fourth order in t. For some, 
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a sixth order expansion was necessary, where we found 
that at least the term B 2 RB 3 helped to distinguish these 
graphs. We tried a sampling of graphs from the families 
up to N — 40, and found that this sixth order term suc- 
ceeded in every instance we tried. 

VII. DISCUSSION 



port from a National Science Foundation Graduate Re- 
search Fellowship. The authors thank the extremely 
helpful HEP, Condor, and CHTC groups at Madison for 
their invaluable assistance in the numerical components 
of this work. They also thank Conrad Sanderson for his 
Armadillo C++ linear algebra library and his help with 
its use, as well as S-.Y. Shiau for his helpful correspon- 
dence on the problem. 



This paper takes several steps towards characteriz- 
ing the additional power that quantum walks of inter- 
acting particles have for distinguishing non-isomorphic 
strongly regular graphs, as compared to quantum ran- 
dom walks of nonintcracting particles. We prove analyt- 
ically that quantum random walks of two noninteracting 
particles, either bosons or fermions, cannot distinguish 
non-isomorphic strongly regular graphs from the same 
family. We investigate numerically the quantum time- 
evolution operator for a quantum random walk with two 
interacting bosons and show that the resulting Green's 
functions can be used to distinguish all non-isomorphic 
pairs of strongly regular graphs that were investigated. 
We perform a much more comprehensive numerical test 
of the interacting particle algorithm than has been done 
previously, and find that quantum random walks of two 
hard-core bosons successfully distinguish all known non- 
isomorphic pairs of SRGs, which include graphs with up 
to 64 vertices, and family sizes as large as 32,548 non- 
isomorphic members. 

We now discuss how our results are relevant to possible 
algorithms for solving GI. If our algorithm for two hard- 
core Bosons does indeed distinguish arbitrary graphs, 
then GI is in P, since the number of particles is fixed, 
and only the lattice size increases with number of ver- 
tices. But, if GI is not in P, then for some pair of graphs, 
our two particle algorithm must break. Hence, if such a 
case was identified, then one could try increasing parti- 
cle number, which could potentially place GI in BQP, the 
complexity class solvable efficiently on a quantum com- 
puter. Unfortunately, we exhausted our test cases (the 
SRGs), so we could not test this hypothesis. 

Although we found that two non-interacting Bosons 
were not helpful for distinguishing SRGs, we suspect that 
at larger numbers of non-interacting particles, QRWs 
might be able to distinguish SRGs. This suspicion is 
due to refs. [34Tl36j . where non-interacting QRWs ex- 
hibit effective interactions in the statistical limit. These 
effective interactions might produce distinguishing power 
comparable to our explicit hard-core interaction, while 
providing for an easier analysis. It would be interesting 
to examine the several-particle non-interacting QRW to 
see if this is indeed the case. 
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Appendix A: Checking the two-particle matrix 
elements 



In this appendix we demonstrate that eqs. 15 and 23 



are equivalent to Eq. [7]for both Bosons and Fermions. To 
show that the Boson matrix elements as given in Eq. [15] 
are correct, we evaluate the three types of basis elements 



we have in eq. 15 When i =/= j and k =/= I, Eq. 15 yields 



s (ij\U 2B \kl) B = 



( + in\' 



V V2 

\kl)'+\lk)' 



-i(I + S)(A©A) 



UK 



V2 
\kl)' + \lk)' 

7T~ 



(A ©A) 



(Al) 



— —Aikbji — Aji5ik — AuSjk — AjkSu. 
When i j but k = I, we find 



,(ij|H 2B \kk) B = 



\ V2 
1 



(I + S) (A © A) + UK 



\kk)' 



V2 



(A © A) \kk)' 



V2 
+ AjkSik) 
2 

~ ~vi 

When i — j and k = I, we have 



(AikSjk + Aj k 5ik + AikSjk 



{A ik 5jk + A jk 5ik) . (A2) 



b (m|H 2 b \kk) B 



-i(I + S) (A©A) + [/R 



U(ii\'K\kk)' 

us lk . 



\kk)' 



(A3) 



These expressions are all exactly the same as we found 
through the definition of the Hamiltonian in Eq. ref- 
bosonelements. 

We can work the same exercise for the Fermion Hamil- 
tonian, Eq. fermionoperatorham as we did for the Boson 
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Hamiltonian in eq. AT The calculation for the fermions 
yields 

F (ij\H 2 F \kl) F = SikAij + SjiAik-SuAjk-SjkAii, (A4) 
which is identical to the result obtained from eq. [7] 



Appendix B: Counting the elements in the 
non-interacting Boson evolution matrix 

To prove that the two-particle walk cannot distinguish 
two non-isomorphic SRGs, we use the method that Shiau 
et al. introduced for one-particle walks; we show that 
all the values and degeneracies GFs, matrix elements of 
the evolution operator, can be expressed of functions of 
the SRG family parameters. For one particle, Shiau et 
al. [3D] considered the on-diagonal and off-diagonal ma- 
trix elements separately. For the two-particle evolution 
operator XJ2B , we perform a similar trick by partitioning 
the matrix elements according to two parameters (a, b) : 
the total number of distinct indices (a) and the number 
of indices shared the left and right sides (b). For exam- 
ple, (34| Ub 1 24) falls into the element class (3,1), since 
it has three distinct indices (2,3,4) and the left and right 
side have one index in common (2). In total, the two- 
particle Boson evolution operator has seven such classes: 
(4,0), (3,0), (3,1), (2,0), (2,1), (2,2), and (1,2), which 
together partition the set of matrix elements. 

Within each of these classes, the various possible ele- 
ment values are listed in table |lj Counting the number 
of occurrences is performed by means of combinatorial 
sums. First, we consider the symmetry class (4,0). Since 
i 7^ j and k ^ I, we use eq. [5]to write the matrix elements 
as 



(ij\U B \kl) = a 2 (6 ik 5 jl +6 il 5 jk )+2j3 2 

+ ^(AikAji +A u A jk ) 

+ a[3(5 ik + S jk + 5u + Sji) 

+ aj(S ik Aji + Sj k Au + SuAjk + 5 iA ik ) 

+ pj(Aji+Au+A jk + A ik ). (Bl) 

One possible value for this matrix element, 4/?7 + 2-f 2 + 
2/3 2 , occurs when Aji — An — Aj k — Ai k — 1. Since 
A is a — 1 matrix, the number of ways this can occur, 

i(4,o)«n is s iven b y 



n {4,0)a = / ) AjjAjlAjkAj 

i<j k<l 



(B2) 



J AjiA u A jk A ik (l - 5 i:i )(l - S M ) 

ijkl 

1 -{j2(A% i -2j2(A 2 ) ij + j:jA ij ), 



where the initial sum is constrained to i < j and fc < I 
since we are working with indistinguishable Bosons, and 
hence have a space of dimension N(N+1) /2. By repeated 
use of eq. [T2j we can use the techniques of the one- 
particle algorithm |30j to evaluate these sums in terms of 
family parameters. The values of those pertinent to our 
present discussion are: 



5> 

ij 

E(^ 3 )*. 

ij 
i 



kN 

N(k - m) + kN (^ - M) + 

N (k 2 + k(n(N + yu — 2) + A 2 - 2A/i + A) 

(iV-lMA-M)) 
kNX 



kN (fi(k-\-l) + k + X 2 ). 



(B3) 



Plugging in these results and using the SRG relationship 
(N-k - l)/i = k(k - A - 1), we find 



«(4,o)a = l/4-N(k 2 (v + l) + k(\ 2 -\(fi + 2) + n-l) 



- 2(N-l)n) 



(B4) 



A second possible value for this matrix element is 3/37 + 
7 2 + 2/3 2 , obtained by setting any one of Aji, An, Aj k , 
or Ai k to zero. As a sum, this means that the number of 
occurrences, ra^ou, is 



1(4,0)6 



= 4^^A // A«A !fc (l- A ik )(l- S ik ) 

i<j k<l 

= AjiAuA jk (l - A ik )(l - S ik )(l - %)(1 - S M ) 

ijkl 

= E^%--E^ 3 )«-E( a4 )" 



= N(i(N-k-l)(k + X-n), 



(B5) 



where the initial factor of four is due to the four possible 
ways to pick the A, and (1 — Ai k ) (1 — 8i k ) constrains both 
i =/= k and Ai k = 0. The remainder of the calculation 
proceeds similarly, and the results are listed in table [I] 



[1] R. Motwani and P. Raghavan, ACM Comput. Surv. 28, [2] R. Aleliunas, R. M. Karp, R. J. Lipton, L. Lovasz, and 
33 (1996). 



12 



C. Rackoff, in SFCS '79: Proceedings of the 20th An- 
nual Symposium on Foundations of Computer Science 
(IEEE Computer Society, Washington, DC, USA, 1979), 
pp. 218-223. 

[3] Z. Trautt, M. Upmanyu, and A. Karma, Science 314, 
632 (2006). 

[4] R. Sessions, M. Oram, and M. Szczelkun, J. Mol. Biol. 

270, 413 (1997). 
[5] L. Kilian and M. P. Taylor, J. Int. Econ. 60, 85 (2003). 
[6] Y. Aharonov, L. Davidovich, and N. Zagury, Phys. Rev. 

A 48, 1687 (1993). 
[7] E. Bach, S. Coppersmith, M. P. Goldschen, R. Joynt, and 

J. Watrous, Journal of Computer and System Sciences 

69, 562 (2004). 
[8] D. Solenov and L. Fedichkin, Phys. Rev. A 73, 012313 

(2006). 

[9] A. Childs, E. Farhi, and S. Gutmann, Quantum Inf. Pro- 
cess. (2002). 

[10] N. Shenvi, J. Kempe, and K. B. Whaley, Phys. Rev. A 
67, 052307 (2003). 

[11] A. Ambainis, International Journal of Quantum Informa- 
tion 1, 507 (2003). 

[12] A. Ambainis, Foundations of Computer Science, Annual 
IEEE Symposium on 0, 22 (2004), ISSN 0272-5428. 

[13] F. Magniez, A. Nayak, J. Roland, and M. Santha, in 
STOC '07: Proceedings of the thirty-ninth annual ACM 
symposium on Theory of computing (ACM, New York, 
NY, USA, 2007), pp. 575-584, ISBN 978-1-59593-631-8. 

[14] V. Potocek, A. Gabris, T. Kiss, and I. Jex, Phys. Rev. A 
79, 012325 (2009). 

[15] D. Reitzner, M. Hillery, E. Feldman, and V. Buzek, Phys. 
Rev. A 79, 012323 (2009). 

[16] G. Agarwal and P. Pathak, Phys. Rev. A 72, 033815 
(2005). 

[17] K. Eckert, J. Mompart, G. Birkl, and M. Lewenstein, 

Phys. Rev. A 72, 012327 (2005). 
[18] C. M. Chandrashekar and R. Laflamme, Phys. Rev. A 

78, 022314 (2008). 
[19] C. Ryan, M. Laforest, J. Boileau, and R. Laflamme, 

Phys. Rev. A 72, 062317 (2005). 
[20] B. Travaglione and G. Milburn, Phys. Rev. A 65, 032310 

(2002). 

[21] H. Schmitz, R. Matjeschk, C. Schneider, J. Glueckert, 
M. Enderlein, T. Huber, and T. Schaetz, Phys. Rev. Lett. 
103, 090504 (2009). 

[22] K. Manouchehri and J. B. Wang, Phys. Rev. A 80, 
060304 (2009). 

[23] A. Childs, Physical Review Letters (2009). 

[24] D. Spielman, Proceedings of the twenty-eighth annual 



ACM . . . (1996). 

[25] P. W. Shor, in SFCS '94: Proceedings of the 35th An- 
nual Symposium on Foundations of Computer Science 
(IEEE Computer Society, Washington, DC, USA, 1994), 
pp. 124-134, ISBN 0-8186-6580-7. 

[26] S. Hallgren, C. Moore, M. Rotteler, A. Russell, and 
P. Sen, in STOC '06: Proceedings of the thirty-eighth an- 
nual ACM symposium on Theory of computing (ACM, 
New York, NY, USA, 2006), pp. 604-617, ISBN 1-59593- 
134-1. 

[27] C. Moore, A. Russell, and P. Sniady, in STOC '07: Pro- 
ceedings of the thirty-ninth annual ACM symposium on 
Theory of computing (ACM, New York, NY, USA, 2007), 
pp. 536-545. 

[28] T. Rudolph, arXiv:quant-ph/0206068vl (2002). 
[29] V. Gudkov and S. Nussinov, arXiv:cond-mat/0209112v2 
(2002). 

[30] S. Shiau, R. Joynt, and S. Coppersmith, Quantum Inf. 
Comput. 5, 492 (2005). 

[31] B. Douglas and J. Wang, Journal of Physics A: Mathe- 
matical and Theoretical 41, 075303 (2008). 

[32] D. Emms, E. Hancock, and S. Severini, Electron. J. 
Comb. (2006). 

[33] D. Emms, S. Severini, R. Wilson, and E. Hancock, Pat- 
tern Recogon. (2009). 

[34] R. Burioni, D. Cassi, I. Meccoli, M. Rasetti, S. Regina, 
P. Sodano, and A. Vezzani, Europhys. Lett. 52, 251 
(2000). 

[35] P. Buonsante, R. Burioni, D. Cassi, and A. Vezzani, 
Physical Review B (2002). 

[36] F. P. Mancini, P. Sodano, and A. Trombettoni, AIP Con- 
ference Proceedings 918, 302 (2007). 

[37] C. Gosil and G. Royle, Algebraic Graph Theory 
(Springer, 2001). 

[38] E. M. Luks, Journal of Computer and System Sciences 
25, 42 (1982). 

[39] E. W. Weisstein, Strongly regular graph 

(2010), URL [http://mathworld. wolfram, com/] 

Str onglyRegularGr aph . html . 
[40] E. Spence, Strongly regular graphs on at most 

64 vertices, URL http://www.maths.gla.ac.uk/~es/ 

srgraphs . html 

[41] B. McKay, Combinatorial data, URL http://cs.anu. 

edu . au/~bdm/data/ graphs . html 
[42] U. Schoning, J. Comp. Syst. Sci. 37, 312 (1988). 
[43] Specifically, it has been shown that graph isomorphism 

were NP-complete, the polynomial hierarchy collapses to 

level two [42]. 



